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The author has pleasure in acknowledging his indebtedness to 
Dr. T. J. Baker for information from which this research originated; to 
Mr. 0. F. Hudson for the preparation of the micro-photographs; and to 
Mr. J. L. Dixon, Bowen Besearch Scholar in the Metallurgical Department 
of the University of Birmingham, for the care with which he has repeated 
and extended many of the earlier observations. 


The o)-Functions, a Class of Normal Functions occurring in 

Statistics . 

By E. Cunningham, M.A., Fellow of St. John’s College, Cambridge, Lecturer 
in Applied Mathematics, University College, London. 

(Communicated by Professor Karl Pearson, F.E.S. Beceived May 29,—Eead 

June 25, 1908.) 

I.— Introductory. 

The present paper originated in an attempt to discover the significance of 
^certain functions developed by Professor Karl Pearson in a memoir entitled 
A Mathematical Theory of Bandom Migration,” and by him called 
o)-functions . They belong to the category of normal functions, and are 

.applied in the memoir named to obtain an expansion to represent a distribu¬ 
tion symmetrical about a point in a plane. The distribution is not fixed but 
depends on a parameter a, the function co 2n being a function of x = r/<x, r 
being the distance from the centre of the distribution. The fundamental 
differential equation is {d 2 /dx 2 ~\- (x-\-x~ l ) djdx-\- 2 (n + 1)} co 2n (x) — 0 . 

In the course of the present investigation it soon appeared that the same 
function led to solutions of the equation of conduction of heat in two 
dimensions for the case of symmetry round the origin, the time t taking the 

place of a 2 . In fact, it was found that if a solution of that equation is 

$ 

sought in the form f(t)<p (r 2 /^), that solution is t~ ^ +1) o^nf 2 /^), n being' 
arbitrary. The function co 2n is equal to e~ r2 l u multiplied by a polynomial, 
and is therefore especially adapted to the solution of the problem of the 
•cooling of an infinite sheet, the temperature at a great distance being always 
zero. 

The next step in the paper is to generalise the co-functions ; and all the 
■solutions of the equation \/ 2 u = du/dt are found which are of the form 
,f(t) cj) ( r 2 jt) ©, where 0 is a function of the angular co-ordinates of the point 
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1908.] a Class of Normal Functions occurring in Statistics . 311 

alone. The result is that, as in other problems, 0 is a spherical harmonic, 
or circular function, in the case of three and two dimensions respectively, 
that f(t) is of the form t ~ n , while </> (r 2 /At) satisfies the generalised equation, 

d 2 y jdx 2 -{-(%+z ~ 1 ) dy/dz + (2n—m 2 z~ 2 ) y = 0, 
where £ = rt~~% or xd 2 yjdx 2 + (x +1 + m) dyjdx + (n + Jm) y = 0, if x = r 2 jAt. 
This equation is solved by the Laplace transformation, and the solution, 
{2 r rri)~ 1 e~ x ^e~ ux (uA- iy + ^ n ~ 1 u^ m ~ n du, taken round a loop contour from infinity 
round u = 0, is denoted by co n ^ m . The asymptotic expansion of co n , m is 
{e~ x (~x) n ~& n /F (n — ^m +1)} 



which terminates if either n—\m or n + \m is a positive integer. 

Thus the solutions of \? 2 u = dujdt in the cases of one, two, and three 
dimensions respectively are 

t-W X ieo n+h j (x), t-^^ m co n , m (x)me, r^x^co ^, Ic+i («) y, c ( 6, <f>), 

yic being a spherical harmonic of order k, which reduce for certain values 
of the constants to simple well-known solutions, viz.:— 

(i) To ( — irt)~h~ x (ioi n = — (ii) to t~ 1 e~ x (for n = m =0), 

(iii) to (for n == -J, k = 0). 

The one-dimensional solution obtained is not new, it practically coincides 
with the function of the parabolic cylinder, or the functions considered by 
Sturm* in connection with the conduction of heat, by Thiele and Charlier 
in reference to probability, and by Hermite.f The two-dimensional functions 
seem to be new, but they are shown to be the equivalent in polar 
co-ordinates of the functions considered by Hermite 

(3 ' /d'x) m (3 jdy) n e~$^ + v z \ 

and the relations connecting the two types of function are obtained. 

It is shown in the paper that any function of x and 6 can be expanded in 
the form 

m r2i r rw 

m COS mO, where (—) 7 ra n% m = fe x x~> ? _ TO cos m6 dx d0 y 

Jo Jo 

and a proof of the convergence and validity of the expansion is given in the 
case where / is independent of 9 } with an indication of the extension to the 
general case. Hence a series is obtained giving a solution of the equation 

* Liouville, I and II. 
t ‘ Comptes Rendus, 5 vol. 58. 
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y 2 u = dujdt , which reduces for £ = 4 to a function given all over the plane,. 


viz. 


m (r 2 /44) cos md, the coefficient a n% m being obtained as ! 
above. This series is shown, by direct expansion, to be equivalent to 
Laplace’s definite integral solution of the cooling of an infinite plate, viz., 


/(£ v) e -^“£> 3 + (r~r,wMt-t 0 ) d%d v /47T 


Various other expansions are given analogous to Neumann's expansions of 
J n (r 2 + p 2 —2rp cos Off in a series of Bessel functions. 

The last part of the paper considers the equivalent functions in Cartesian 
co-ordinates, in terms of which Hermite has expanded any function., 
It is shown that an equivalent method of obtaining Hermite’s expansion 
is that of equating the successive moments of the series to those of the given 
function, with special reference to obtaining the best fitting surface for 
a given frequency distribution in two dimensions. The series of approxima¬ 
tions so obtained is compared with that of Edgeworth,* which is ultimately 
the same, and is shown to converge rapidly under the same conditions. An 
example is then given of the contour lines of a surface obtained by taking, 
terms up to the fourth order. These contour lines show much similarity 
with those actually observed in statistical work, such as those given by 
Perozzo in his memoir.f The possibility of obtaining a fairly close approxi¬ 
mation to a given distribution by means of a knowledge of the product, 
moments up to the fourth order is thus established. 


II.—- The Derivation of the Generalised wffunction. 


1. The Fundamental Differential Equation. —In order to obtain the types- 
of function suitable to the problems mentioned, the most general solution 
of the equation y 2 u = dujdt will be found in the form T.K.©, where. 
T is a function of t alone, E a function of r 2 jt , and © of the angular 
co-ordinates. To begin with, the cases of one, two, and three dimensions are 
considered separately. 

(a) One Dimension.— Putting z = r 2 jt in the equation, it becomes 



If u = Z. T, Z and T being functions of z and t respectively, we have 



( 1 ) 


so that both sides must be independent of t and z. 


* 4 Camb. Phil. Trans.,’ August, 1904. 
t 4 Annali di Statistical 3 a, 5, 1883. 
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Hence t~+nT = 0 and 4z~ +(z+2)^+nZ = 0, 
dt oz 2 oz 

where n is any constant. Thus T = t~ n . 

( [b ) Two Dimensions .—The equation in polar co-ordinates is 

d 2 u 1 da 1 d 2 u __ da 
dr 2 r dr r 2 d6 2 dt 

If u = ©y satisfies this equation, © being a function of 6 only, and v being 
independent of 6 , the functions © and v satisfy 



m being' a constant. Thus © = A sin md + B cos m6) and again, putting 
z = r 2 jt , and v = T. Z, the same procedure as above gives T = t~ n and 



( 2 ) 


(c) Three Dimensions .—Exactly similar procedure in this case leads to the 
conclusion that, if u = T. Z. ©, © being a function of 6 and <£, © must be a 
spherical harmonic—say of order k, T must be t~~ n , and that then Z satisfies 
the equation 



(< d ) Three Dimensions—Cylindrical Co-ordinates (p, 6 , £).—In this system 




equation 



The equations (1) to (4) can be brought under the same form. Putting 
Z = 3*X, # m X, 2 * /c X, a£ m X in the four cases respectively, and also z = 4a?, 
the equations become, using dashes to denote differential coefficients, 


xX" + (x + $) X' + (n + $)X = 0 , 
xX" + (x+l + m)X' + (n+lm)X = 0 , 
a>X" + (a? + |. + *) X ; + (» + i*)X = 0, 

^X // + (^ + l + ^) X' ■+ (n + k 2 + lm) X = 0. 


(5) 

( 6 ) 
CO 
( 8 ) 


Taking (6) as the standard form, the particular solution dealt with in this 
paper will be denoted by co n „ hm : the corresponding solutions of (5), (7), and 
(8) are therefore 
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and the complete solutions of the original equations are 

(9) r<»+%Wn, i (*), (io) r<»+%4™ ™ me m (*), 

(11) r(”«)# B8 . u+S («) jfc (0), (12) (a) fe ^ m0. 

The simplest solutions are given in the four cases as follows:— 

(i) n = —(ii) % = 0, m = 0, (iii) w = J, ft = 0, (iv) n = m — ft = 0. 

2 . The Solution of the Fundamental Equation, —-Putting X = e~ vx \/dv 

J c 

and forming the Laplace transformation of the equation, we find that Y must 
satisfy the equation 

v 0-1) dY/dv = {(m—1) v —0 + |m—1)} V, 
giving Y = Av n+ i m ~ 1 (v—l)% m '~ n r 

while the contour c must be such that Ye~ vx (v—v 2 ) has the same value at the 
two extremities. Assuming x to be real and positive, this will be attained by 
taking a contour consisting of a line following the real axis from v =s oo to 
v = a(a<l), then a circle about v = 1 of radius a , and a straight line 
returning along the real axis to infinity. A similar contour encircling the 
point 0 instead of 1 will be denoted by 7 . Putting u = v — l and A == 

we have for co, hm the expression (^7ri)e~ x i e~ ux {u + l) n+ i m ui m ~ n ~hkL 

Jyo 


The integral gives the asymptotic expansion for a> nt m 

e~ x (—x) n ^ m P | (^ + jm) (^—j-m) 

^ ^ L % 

(% + ^m)(% + Jm — 1) {%—\m)(n— 1) 

+ 2 ! .(T 2 


}■ 


which terminates if either n—\m or n + \m is a positive integer 

The most important case is that in which n — \m is a positive integer; in 
that case 

(14) co n , m = e *T(»+ im+ 1){t ( to+ + 2)1!(m-Jm-1)! 

+ ... f — )»“*« - - -. - X . I 

r (n 4- -Jm +1) (A—Jm) ! j * 

In particular, if m = 0 this reduces to Pearson's function aj 2 n* 

When n—\m is half an odd integer, becomes purely imaginary. This 
case occurs below, but no difBculty arises on that account. 

Especially important is the case in which both n~\m and n + \mT are 
integers, for then we have the relation, 

(15) (— (n+ Jm + 1 ) = (— x)~^ l €o n ^ m /T(n— + 
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1908.] a Class of Normal Functions occurring in Statistics. 315 


3. Relations connecting Successive Functions .— 

(a)Difference equations— 

27 tl x dco n _ i ? m /cfe = — f (u + ^ g u 

J y 

= — [ 0/0?6 {^+ clu 

Jy 

— 27T2 ( (^ -f~ ~2^F) CO n } m (^2 2 ^^) i ? }}£ j"• 

Thus x dco n - 1 5 TO / dx -f (72 —|-m) 1 , ™ + |m) co n , m . (16) 

Applying this twice and simplifying by means of the fundamental 
differential equation, we find 

(? 2 -f-|m) o) n ^ m ( 2 t 2 x 1) w n - ljm +(22 — Jm — 1) w„- 2)m = 0. (17> 


(6) The Expression of o) n ^ n as a Differential Coefficient .—Since the 
derivative of any solution of V 2 W = 3 m/ 0£ with respect to t is also 
a solution, and since 3 fit {t n f(r 2 ft)} = t n ~ l $(r 2 ft), we deduce at once,, 
after adjusting the numerical coefficients, the following relations in the casein 
of one, two, and three dimensions respectively (vide (9), (10), (11))— 

dfit {t~ n x%(D n -$^ (x)} = —nt~ n ~ 1 x^co n+ i^(x), (18) 

d/dt {r n ^ m w n _ 1>m («)} = — (n — \m) t~ n ~ 1 !^ m (o ntm (x), (19) 

d/dt {t- n w k w n _ lk+i (i »)} = -{n-\h)r n ~'xl k co n _ h , c+h (x). (20) 


III .—The Linear Function co n+ ^ $. 

1. Restricting the work in this section to the functions arising in the 
one-dimensional case, which are briefly called the linear co-functions , it will 
be convenient to write co v for a> n -i,b and <*'v for -j,. and v in place of n* 
From (18) we then have 

t- v -h*co v = (-B/dtyit-teaoix)} r(i)/r(i; + i) 

= (-y~Hd/dty(t-h - x )/r (v + i) = (-y~ h (0/0r)^ (t-*e-+t% 
assuming that is an integer. If j/ —J is an integer, we have, similarly, 

t- v ~wco v = (-y-^d/dty-^it^x-^ix^/riv+i) 

= ( — y-% (d/dr) 2v ~ 1 {^t~%re~~ r<i l u } = ( — ) v ~*(djdr) 2v (t~h~~ r2 l u ). 

Thus in either case 

= (-)■'-**■' (0/0r) v (e-^)/rV {J (i/ +1)}. (21), 

Thus, as a function of r, wi r is identified with the function considered by 
Hermite.* 


* ‘ Comptes Eendus, 5 vol. 68. 
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2. An Expansion in a Series of Linear a>-Functions .—The expansion of 
ai p, e- x (-xfl<r-v f 1 v(y — V) . v (v — 1) (v — 2) (v — 3) 1 , - f • 

ry^Fiyi 1 ~4T~ + 2!(4xf whlch ’ 11 * 1 

.an odd integer, becomes 


is 


! 


4x 


2'-T{J( v + l)} liO'-l) 1 i (V—3) ! 3 ! 
.and if v is an even integer, 

- x v k-®)-* r l Ax 


+ ... 




r I 


‘FT{\(v+l)} 


If, 


i ( v —2)! 2 ! 


+ ...- 


(4o?)* v *) 
“j* 


Hence, = coefficient of in v ! e^ VMq -^ x ~ x j{—xf2 v Y {^(z/ + l)}, 
and therefore 

.»-**■-* cos 2qs/x = it(2qyr{i(y + l)} co hv /T(v + l) = ctq^ v T(i)/V (±v + 1) 


and 


( 22 ) 


.arM* x sin 2q\fx = 2(2y)T{J (v + 1)} oo iv /T(v + l) = Xq v co iv T(^)/r(^v + l), 

(23) 

the former summation being for even values of v and the latter for odd, or 
together, 

x~h^^ = \Z~^ri(q) v co hv / r(Ji; + l). (24) 

v — 0 


3. Application of the above Expansion .—-The foregoing properties will now 
he applied to obtain an expansion of Laplace's solution of the problem of the 
^cooling of an infinite bar. We have 

. {2?r(^-4)}-^“ (r -p )2/4(f -^ = (27T)" 1 I* e-vW-to)~^'~p)du 

J -GO 

= (An) ~i i f e ~ uH - iur p<fl u » tn v ty~V><o iv (p 2 /At 0 ) du/V (\v+ 1), 

J — 00 

= (A-rr)-ip^X{-t 0 )^-»a) iv /T^v + l)(0/3r)-f e~ u2t ~^du, 

J —CO 

= (2f 0 )-* (p»/4«6)* (r»/40* $ ai, (p 2 /4< 0 ) o>„ (r 2 /40 

r{i(r+i)}/r(^+i), 

= ( — 2t 0 yht %ixi % (t 0 Jt) v+ ia)' iv (£) &>I„ (*), 
using (15) and putting £ for p 2 /At 0 . 

r® 

Now, as Laplace as shown, Y = /(£) (47 rtyh-^-^d^ is a function of cs 

J -GO 

satisfying 0 2 V/3^ 2 = 0V jdt and reducing when t = 0 to fix). Putting (£—£ 0 ) 
in place of £, and expanding the subject of integration as above, we obtain 
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the following series to represent a solution ['of the same equation reducing 
when t = t 0 , to fix) • 

%Ov (tolty^Ww (x), (25) 

where ci v — —f f f(p)ef-co' v (^)d^. 

J —oo 

In accordance with a remark made above, it may be noted that, when ay v is 
real, a>' v is imaginary, so that the expansion is wholly real. A different 
method of establishing the same expansion is given below. 


IV .—The General ay-Function—Integral Properties . 
Let 

f 00 

e x co ni m (x) &v s _ m (x) dx 
o 


Jo Jy Jy 

The asymptotic expansion of the ^-function shows that the integration in 
regard to x to the infinite upper limit is valid, and the order of integration 
may be inverted. Thus 


— 47t 2 I^ n t }m = dul dv(u + l) n+ * i ^- n - 1 (v+l) nf -* , ty-i* l - n, - 1 (u + v+l)r\ 

Jy Jy 

For every value of v the subject of integration is zero to the second order 
at least when u is infinite if m is positive, and vice versd if m is negative. 
Taking m positive for the sake of argument, 

47T“I w? m 

= — f du f dv (u + 1 yt+hm u lm~n~\ ^ ^ _|_ V + 1)“L 

j —j/—1 J y 

= ( — ) m 2m f dv v n ~ n '- 1 (l + v) n '- n - 1 = (—)»( —47t 2 ), 

Jy 

if n = nl\ or = 0 if n + n\ Thus 


but 


-co 

j e x (Oyi t m (<r) g)h> t — m (pc) dx — 0, 

Jo 

/»CO 

I S x CO n ^ m (pc) G) n: — m (x) dx = ( l) m . 

Jo 


(26) 

(27) 


Again, consider the moment integral 


27riM r = 


= 27 rii (v n , m ( x ) x r dx = f x r dx\e ( u+1 ) x (u + l) n+ % m u& n ~ n ~ 1 du 

Jo Jo J v 

= r (r+l) f (m+1 ) n+ i m ~ r - 1 ui m - n ~ 1 du. 

Jy 
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Thus, if is a positive integer, (l + w) w+ ^ m ~' r ~ 1 is a polynomial and 

the subject of integration has no residue at n = 0. But if n + ^m—r is zero or 
a negative integer, the residue is (—) n -+1 — n — \m )... (r—m)/ (71 — |m)!. 
Thus 

M r = 0, (—(—)»“*» (~-) M “^r(r+l)r(r— m+l)/(n—%m)\ 

T(r-n-±m + l), (28) 

according as r <, =,or > n + ^m. In this r is not necessarily an integer, 
but n—\m and n + ^m — r are assumed to be integers. 

If these results, (26) and (27), be applied in the one-dimensional case 
considered above, assuming that /(£) = we obtain at once 

a v = — I Consequently, if £ = r 2 /4£ 0 and x = r 2 /4:t, 

Jo 

%a v (t 0 /t) n+ ^co v (x) is a function satisfying d 2 u/dr 2 — du/dl, and reducing 
when £ = to to the given function/^). Thus the expansion (25) is obtained. 

Y .—The ay-Functions in a Plane . 

Coming to the special values of n and m suitable to the development of 
functions in a plane, m is an integer and n—\m is an integer. Using 
relation (19), we have 

On, m = t' (n+1) ( r 2 /t ) -**©«, m (r»/4^ ^ 

= ( — d/dt) n {t~ l (r 2 /t)~% m a)Q tm e ime } T (^m+ 1)/T +1) 

= { — y 2 }^{r 1 (r 2 /t)~* m G)Q im e Lme } r (Jm + D/r (Jm + w+1). (29) 

In particular, if we start from the function independent of 0, i.e. put m = 0 , 
O n?0 = r n+l co n (r 2 /4:t) = (-3/30” (rh-^/nl, so that 

= (—<) n /! (3 2 /3r 2 +1 jr 3/Sr)^ - ^ 4 *. (30) 

This particular set of functions may conveniently be called the zonal 
^-functions. 

If m be not zero, but n—\m a positive integer, so that a n > m is polynomial, 

. n*,» = {-d/dty-i™ {r«« +1 > (?/t)-**a im , m e me } 

= (—d/dt) n -l m { r m t -m-i-ryit e an8} r(m+l)/r (|m+%+1) 

= ( — )^“s m (y 2 )n-|m ^ r m^-m-l e -r2/4^im0j. p (^-f 1)/T (|m + % ■+ 1). (31) 

Thus all the functions found are expressed in terms of the derivatives of 
the fundamental function with respect to r, 6 , and t. 

Since the derivatives of any solution of v 3 ^ = du/dt with respect to 
Cartesian co-ordinates x and y are also solutions, we may expect all solutions 
to be expressible in terms of the derivatives with respect to x and y of the 
fundamental solution. 
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The following equations may be shown to hold for the planar w-functions, 


= ( B/Bx + iB/By) m (t 


and, therefore, n—\m being an integer, 


Q n , m = ( 3 2 / Bcc 2 + B 2 /3y 2 ) n lim ( B/Bx + iB / By) m ( t x e ~ r2/4i ) 
= (3 /Bt) n ~i m (B/Bx+iB/B y y i (rh-^% 


or (3/ Bx + id / By) 71+ (8 /3x—i3/By) 71 (t 1 e, r2/4 *); (32) 

e.g., we have 


{B 2 jBx 2 —B 2 jBy 2 ) (t-h-'W) = (d^ 2 )- 1 (r 2 / 1) e~ r ^ 4t cos 20. 
2 3 2 /3x3y (t~ 1 e~ r2 ^ t ) = (d^ 2 )” 1 ( r 2 / 1) sin 20, 


and (3 2 /3x? + 8P/3y 2 ) (t 1 e~ r ^ u ) = — (t 2 ) 1 (1— r 2 /4:t) e" r<2 / u , 


giving the values of B 2 /Bx 2 , 3 2 /3y 2 , B 2 /Bz 2 in terms of fli 0 and Oi 2 . Thus the 
connection is established between the planar w-functions and Hermite’s 
functions B r+s /Bx r By s ^e~ a ^ x2+ ^}, to which reference is made below. 

VI.— The Expansion of any Function in Planar co- Functions. 

1. Suppose that a function of r and 0 may be expanded in the form 
/ (r, 0) tta nm (xf m (o nm (x) cos m0 + ttb n , m (x) sin m0, 


n m 


where x — r 2 /4ct, t being a parameter. Then, using the Fourier method for all 
values of m, we have 



and 


f(r, 0) sin m0 eld — nrXb ntW x^ m eo niVn . 


o 


In order to find the coefficients a n%m% b n , m% we now use the integrals 
developed above, and obtain 



with a similar expression for l n%m . Thus, in order to obtain the coefficient 
of we make use of the allied function o) n ^ m , which bears to it the 

relation (15). 

The values of m and n required to be taken do not appear, but it will be 
shown presently that m takes all integer values and n takes those values for 
which n—\m is zero or a positive integer. The expansion is therefore one 
in a series of polynomials. The corresponding expansion in Cartesian 
co-ordinates is given by Hermite.* 


Y 2 


* ‘ Comptes Bendus,’ vol. 68. 








Downloaded from rspa.royalsocietypublishing.org 


320 


Mr. E. Cunningham. The Functions, [May 29, 


2. The Expansion of the Laplace Definite Integral in a Plane .—An expan¬ 
sion analogous to (26) will first be obtained, x standing, as everywhere, for 
r 2 /4 t, let z stand for re te j 2 ft and z for re~ l6 /2\/t. Then, taking all values 
of n and m for which n—\m and n + \m are both positive integers, 


n ,m = e x (n + %m )! -f 1 
L m!(n—i 


( n—rn )! (rn +1)! 1! (n— —1)! 


+ • 


■Am)! (n—lm) ! J 


(n + -|m)! (n —|m) , 

= (n + \m )! x~i m . coefficient of h n ~^ m k n+ ^ m in e hk+ ( k ~ h ) ^ x ~ x . 

oo m — +2 n 

or t t m h n ~i m k n+ i m /(n+ Jm)! = e< A+ < 4- ^ x \ (34) 

n = 0 m = — 2 n 

% taking all half-integer values from 0 to infinity. 

Again, taking the same range of values for m and n, 

^ x i m e mi % nm h n ~i m k n+ ^ m /(n^^m) ! = Vx) __ e (h+z)(k-z') t ( 35 ) 

Expanding the last expression by Taylor’s theorem, we have 

x* m e mid co nm = ( — )»+*»» (3 /dz)*-*™ (3/3 zy+* m {e~ zz ')l(n-\m) !. 

Thus, putting £ = re 10 and 7 ] = re lQ \ 
x^ m e mie a) nm == (-f+Hr +1 (d/d!;) n -& n (d/dn) n+ % m (rie-*/")/(n-|m) L 
This gives the proof of the relations given above (32). 

Putting th = p-\-iq and —rfc = p—iq in (35), that expression becomes 
e p^+ g ^+i(pu+qv)-x^ w here (uv) are Cartesian co-ordinates. 

Hence 

f oo /» 00 

dp dq 

—00 •* — 00 

= (27 t) _1 f f °°<2p clqe-r +q2)t - lipu+,lv) h n ~i m k n+ i m /(n + |m) f 

J —00 J —00 

by (35), where % = (% 2 + -ri 2 )/4io. Now put A' = h</(t/t 0 ) and /</= k^/(t — to) 
and A' =y> + t2', k' = p—i,q. Then the above becomes 

22 (_)»-i»* (3/3’) n ~i m (3 /dz') n+ i m 

f oo /.oo 

g-(p 2 +s 2 ) t-*(px+q!/) c ip c iq t 0 n t- n ~ l e^Xii m e m ^co nm (*i)/(w + -Jm) ! 

— 00 J —00 

= tt ( —<%, m (a?i) t 0 n r n -\ 

The summation throughout is for all values of rn and n such that n 4- \m 
and n—\m are positive integers. There will thus correspond to each term in 
the summation another with an equal but opposite value of m, and, using (17) 
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these terms will only differ by being changed to e~ mt ^~ e K Remem¬ 

bering also that n—\m is an integer, ( — )~ 2n = (—)~™. Hence 
( t—1 0 )~ h~ (*- w 

oo 

== 2 C Xl G) n t o (x') CO'fi t n ^ 

n = 0 

+ 222 (—) m e x ^ m xi~^ m a) nm (x) (o?i) cos m ($--0) to n t~ h ~\ (36) 


This expansion is suggestive of Neumann’s expansion in the theory of the 
Bessel function. Another theorem is given below, which, in some respects, is 
a nearer analogue. The development of the Laplace definite integral follows 
at once from (36) and is identical with that in the last paragraph (33). The 
method of this paragraph shows that the values of m and n are as there 
stated. 

3. The Expansion of G> w ,^{(r 2 —2rp cos (<£ — 0) + p 2 )}/4:t in . a Series of 
ay-Functions .—As in the last section, 


1*00 -00 

= (27t) -1 dp 

J —00 J — 00 



dp dq 

ZZ (— h) n -i m k n H»y»-!“£>•+i>»/( K _| m )t (% + lm)! 4“ 


(the summation being taken for the same values of m and n as above) 


= (2 7 r)- 1 f CO [” dpdq(d/d%y-l m (d/dri) n+ i m e 

J _oo J —oo 


- (^2 -f gr2) t— t (p U + g V y 

p2n e mi$ j(ji— |-m)! (n + Jm)! 


= 22 (—) w *H m 4 H n m e mi (<f>~Q)p 2n / (%—-|m)! 

= 22 (—) w+ ^~ 1 (p 2 J4t)~ n x~2 m (o n ,m(%)cosm(<fi~0)/(n+ini)l 
where the summation extends only to values of n such that n—\m is a 
positive integer. Thus G)o ;0 {(r 2 —2rp cos(<f>—0) + p 2 )/4:t} is expanded in a 
series of co-functions in the variable (r 2 /4t\ and from this maybe derived, by 
means of relations (32), the expansions of co n , m {(r 2 — 2?y>cos(</> — 0) + p 2 )/4:t} } 
using differentiations with regard to t } x h and y\. 

4. The Convergence and Validity of the Expansion of a Function in a Series 
of a-functions .—The proof of the validity of expansions in a series of normal 
functions is always a difficult matter, and a general theorem to cover them 
has yet to be established. A discussion is here given of the expansion in a 
series of zonal co-functions of a function symmetrical round a point, which 
seems capable of extension to the general co-functions in a plane without 
much increase of difficulty, though with considerable addition to the analysis. 
The multiplication of different discussions for each particular set of normal 
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functions seems undesirable, however, and the present investigation will 
therefore be limited to the single case of the zonal functions. A different 
method appears possible in the discussion of the validity of the processes in 
the development of the Laplace integral given above. 

In the following, the basis of the work is the definite integral expression 



for the ^-function, viz., {2m) 1 e(u + l) n io n ~ l du . 


y 


The expansion to be justified is 


f{x) = ?La r co r (x), where a r — | f{y) w r {y) e y dy. 



Consider the sum of (n+ 1) terms of the series 


ta r 0 ) r (x) = S I f(y) e y ( 0 r (%) a> r (y) dy 


0 0 J 0 



0 



{^ + l) n+1 (v + l) n+ hr n ~ 1 v- n ~ 1 --l} e-w-™ = -(47T 2 )” 1 {In+1-Io}. 


A necessary condition is at once obvious for the validity of the processes 
involved. For, if y is very large, since the ^-contour includes points for which 
the real part of u is negative, e~ uy is not everywhere finite on that contour. 
But if / {y) be such that, no matter how great y may be, / {y) < e hj , where k 
is some finite negative constant, we may take the ^-contour such that the 
real part of n is always greater than k, and then the integrations required are 
justified. 


Let 



J r (x,y)= e~ u y~ vx (ii + l) r (v + iy u~ r v~ r (tt + v + l)- 1 diodv. 


Then, if y is finite, J 0 is zero, since the ^-contour includes no point at which 

vanishes for any value of v on the contour 7 . Also, by the 

condition assumed to be satisfied by f(y) 9 the part of I 0 arising from values of 
y greater than a finite magnitude r] can be made as small as we please by 

taking rj sufficiently large. Thus To = 0. In treating I n +1 it has to be 

remembered that the limit of the integral is sought when n tends to infinity, 
and that since the contour 7 encircles the point 0 , and the two contours for 
u and v have been assumed such that u + v+l does not vanish for a pair of 
values within these contours, | {u-\-l)ju | and | (v + l)/v | are both greater 








Downloaded from rspa.royalsocietypublishing.org 


1908.] a Class of Normal Functions occurring in Statistics . 323 


than unity on the major portions of the contours, so that the subject of 
integration therefore becomes indefinitely great as n increases. 

r 00 

Consider now the integral J r (x, y) cly and divide into two portions for 

Jo 

which y is greater and less than x respectively, calling them Ki and K 2 , 

Let K a be the same integral with the limits a and x (ct < x ). 

Then 


K a = j dy e x f du j 

J a J v J * 


dv e~ m ~ w J (u + l) r (v+ l) r u~ r v~ r {u+v+ 1) -1 


= — j du^ dv e~ x( ' 1+u+ ^ (u+l) r (v+l) r u- r ~ 1 v~ r (u+v+iy 1 

+ ( du | dve~ x(1+v )- m (u+iy(v+iyu-’^v-r(u+v+iy 1 
= — iK a + 2 K a (say). 

Consider the latter part ; it is absolutely and uniformly convergent, 
and may therefore be written 

—djdxi j* du dv e~ x ^ v ^~ a \ii + V) r {v -f l) r ~ht~~ r+1 v~ r (u '+v + 1)*"' 1 = —3L r /3$ 
J Y J Y 

(say). Within the contours considered there is no pair of points for which 
■(u + 1) = 0. If the ^-contour be extended so as to include all the points at 

which (u + v + 1) vanishes for points on the t'-contour, we add to L r the integral 

f f e~ x ( 1+v )~ au (t{,+ iyu~ r ~ l (v-i-iy* 1 v~ r (it+v+l)* 1 dti dv, the contour 7 , 

J yJ y] 

including u = — v— 1 , but not u — 0 , and this addition 


= -2ttH e ~< x ~ a ^ 1+v \l +v)~ 2 dv = 0, 

Jy 

since the contour 7 does not include the point ( — 1 ), and since x—a is 
positive. Thus the ^-contour in L r may be supposed extended so as to 
include the point at which (u-\- v-\- 1 ) vanishes for all values of v or 7 . 

Let the v-contour be now specified more precisely as that part of a circle of 
large radius E, with centre at v = —that lies on the positive side of the 
straight line on which the real part of v = — together with that part of this 
line which lies within that circle. This contour is considered in three parts, 
a, 7 , as follows: a, the straight portion; / 3 , that part of the curved portion 
adjacent to a and subtending a small angle 8 at the centre of the circle; 
7 , the remainder. 

For values of v on a, the ^-contour, extended to include — (l + -y), is taken 
to be a contour just surrounding the complete ^-contour. For values of v 
on /3 and 7 this contour is extended in the negative direction so as always 
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to be a short distance on the negative side of u = —^—1. The ^-contour 
is considered in three portions similar to a, /3, y, viz.: a h the straight portion ; 
j3i, the arcs contiguous to and subtending an angle 8 at the middle point 
of ai ; 71 , the remainder. Thus the whole of the double integral is divided 
into nine portions obtained by combining one of the partial contours a, / 3 , 7 
with one of the partial contours a h fii, 71. 

Now e~ r!v (\. + v~ l ) r remains always finite as r and v become large, except 
for v = —1, a value which lies outside the region under consideration. Thus 
je“ aw (l+*T 1 ) r | < A\e~ xv+r / v \ and le~ au (l + u~ l ) r \ < B \ e ~ au+r l u \. Thus, if 
v = He LQ \e~ xv (l + v^yi < A^ -2:R+r/R)cos ^, and, if u = RV^, \er au (1 +%C l ) r | 
< Be(-« R '+»*/ E, ) co s^ # Hence, if 

r/R ' 2 = c < a < x, \e~ au (l+u'~ 1 ) r \ < Be~^- C > R ' C0S *, 
and therefore, on the arc 7 , cos </> being > S } 

\e~ ou (l + 'ZT* 1 ) r | < Be~ (a ~ c ) R ' s < Be~< a ~ c ) pS ~* (where R' = 

and this quantity tends to zero as 8 diminishes, provided p remains finite. 

On the arc fii, the expression certainly remains finite. 

On the straight part a h | (u+l)/u \ takes its greatest value at the point 
where it adjoins i.e. {(u + l)ju} r is finite on u lt but along a length (&R')£ of 
a measured from the real axis it becomes comparable with (1 — l/&R') r , i.e. 
with (1 — l/yfcR') cR/2 , i.e. with e~ cR \ where R' is large and c finite. Similar 
considerations apply to the various portions of the ^-contour, and putting 
them together we deduce that each of the portions of the double integral 
tends to zero as r tends to infinity. Thus l 4 = 00 L r = 0 for all values of x, 
and therefore = —dL r /d x also tends to zero as r becomes infinite, for any 
value of a less than x, and in particular if a = 0 . 


Thus 



e -x(i+v)-uy (u + iy (v+iy u" r v~ r (u + v + YyHu dv dy 


e -(i+u+v)X + ly + iy u -r-i u -r (u + v+ iyi du dv , 


and is independent of the lower limit a. 

The u and 'y-contours being the same, this last integral will be unaltered if 
u and v be interchanged. Hence 

K a = e~ x ( 1+w+t, )(^+ l) r (v -f l) r u~ r v~ r +1)" -1 (^ —1 -f- du dv 


— — 1 . e -z(l+u+v) (^_|_l)r (v + iyu- r ~ 1 V-‘ r - 1 { 1 — (u + v + 1) 1 }dudv 



y J y\ 


= i {H+9H/3*}, 
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where 

H = f f e x ^ l + u+v ) + (y+\y u r h) r 1 (m + v+I) 1 dudv 

JyJyj 

= - I e~ x ^ 1+u+v ^u^l) r (v+l)Ho~ r ~ 1 v~ r ~ 1 (u + v + iy'dudv + e, 

JyJ y x 


where € tends to zero as r becomes large ; since, by reasoning similar to that 
used above, the double integral H with the ^-contour extended to include the 
point —(1+u) vanishes as r tends to infinity. . 

Thus 

Lt r = O0 H = — 27m — u~ l (u + l)~ l du = — 47 r 2 and U r - ^dH/dx = 0. 

i y 

Finally, therefore, Ii r = 00 K a = — 27 t 2 , whatever the value of a . 

Apply now the same treatment to the integral K bf viz.:— 

rb 


m e *( 1+w ) , u y (u-\-l) r (v+Vf u r v r (u + v + 1 ) dudvdy 

y 

r= f f e~ x ( l+u+v ) (u+l) r (v + l) r ir r ~ l v~ r (it + v+l)" 1 du dv 

JyJy 

— f I e~ x ( 1+v )- bu (u+iy (v + l) r u- r ~ 1 v- r (u + v+iyWudv. 

JyJy 

The coefficient of u in the exponential in the second integral being now 
greater than that of v, we carry out first the integration in regard to u. The 
first part of K$ is, as above, equal to 27T 2 . Thus 


K b = 2tt 2 — f f e- x ( 1+v )- bu (u + iy(v+iyu- r - 1 v~ r (ii + v + iy 

JyJy 


l du dv. 


In extending the 'y-contour of this last integral, instead of the ^-contour 
as in K fl , we add to it 

I du e ~ x ( 1 +v)~bu l)r ty+iy (U + V+1)~ 1 dudv 

Jy Jy x 

= 2irt f du e~( b ~ x ) u /u = —47r 2 ; 

Jy 

71 , as above, standing for a contour encircling — (1 + u), and not zero. As 
before, the extended integral is zero, and therefore the unextended integral 
equals 47 t 2 Thus K* = 27 t 2 —47t 2 = —27r 2 for all values of b. 

We have 

Lt r = oo { X f(y)dy f f e ~ u ^~ x(l+v ) (u+l) r (v+l) r %r r v~ r (u + v+l'y 1 dudv 

JO J yJ y 

= -2 nr f(x~) 
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and 

Lt r = „ j" f(y) cly j j e - w J- x ^ v ^ (w+1) 1 ' {v + l) r u~ r o~ r (u+v+ l) -1 du dv 

= 2^f(x+)y 

and, therefore, finally 

n 

Ltn^nSdrCCrX = |{/ ( X ~) +f( X + )}> 

0 

and if the function / is continuous}, 

n 

(#y) —• f (^)* 

0 

It will be noticed that the success of the above method depended on the 
subject of integration in the beginning being a geometrical progression. 
If the general expansion for the unsymmetrical case (33) is treated similarly, 
the subject of integration is a double geometrical progression and may be 
treated similarly, but the analysis will not be carried out here. 

VII.— The Application of the co-functions in the Theory of Probability. 

1 . As was mentioned at the outset, the linear ^-functions have been 
recognised for some time as suitable functions for the purpose of representing 
a given frequency distribution to any desired degree of approximation. 
The method commonly adopted in fitting the coefficients is to make the 
successive moments of the required series equal to those of the given 
distribution.* But it does not appear to have been noticed anywhere that 
this method is exactly equivalent to Hermite’s method of obtaining the 
coefficients. This is an immediate consequence of the fact that the function 
$ x% u n is a polynomial. Thus the practical method by moments, which 
commonly appears as an artificial way of obtaining a series to represent a 
distribution of a particular type, is really the exact and unique way of 
representing an arbitrary function by a series of co-functions. Edgeworthf 
has another method of obtaining a series of successive approximations to a 

r—t 

given distribution in the form exp {%k r { — d/dx) r jr !} {(27 rk)~^e~' x ^ 21 '}, which, 

r — s 

if the formal operator be developed, leads to a series of w-functions. The 
method of obtaining the coefficients k r is that of identifying the moments, so 
that ultimately the series as obtained by Edgeworth must agree with the 
direct expansion in co-functions. It might be, however, that, for statistical 
distributions, the method of development used by Edgeworth would give a 

* Charlier, loc. cit. 

+ ‘ Camb. Phil. Trans., 5 Aug., 1904. 
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more rapid approximation. The argument as to the magnitudes of the 
coefficients k r> however, appears equally applicable to the present series.: 
It is shown in the memoir quoted that the coefficients k r rapidly diminish 
only if the standard deviation be small. Let this be assumed in forming 
the series a H u n (x) f where u n — (0/0^) w {( 27 rcr 2 )~ie~* 2 / 2<72 }. Then, if /jl h be the 

nth. moment of the distribution, and ir r = e~ u ^ho r du y we have 

J — 00 

n 

fin = 2 ( — Y Cl r <T n ~ r 7T n— jOfl ! /(n*-r)l (2ir)^ 

0 

and in particular fi 2 = (27r)~* {aoa 2 7T2 + 2<x 2 7r 0 }, 

and hence, since <r is the standard deviation, nv 2 is zero. How, following 
Edgeworth’s argument, fi n is of the order of magnitude a n , and hence, by 
taking the above equations for different values of n, we see that the 
coefficient a n is of the order a n —the word order being used in the same 
sense as in the work cited. Thus the coefficients a n will diminish with about 
the same rapidity as the quantities k n , while the formation of the equations* 
is much simpler. 

Pass now to the case of two dimensions. Reference has been made to* 
Hermite’s note on the expansion of a function in a series of the form 
tta mn d m+n /dx m dy n e“* ^ 2+2/ ^ + W, 

which has been shown to be equivalent to an expansion in co-functions. It 
is a problem of some importance to obtain a series to represent a frequency 
distribution in two correlated characteristics, such as, for example, is con¬ 
sidered by Perozzo in his memoir on the ages of husband and wife at 
marriage.* The lines of equal frequency given by Perozzo are far from being 
concentric ellipses, as they would be if the frequency surface were a Gaussian 
normal surface z = hx y+W\ It is therefore worth enquiring whether 

a few terms of the expansion in co-functions can be made to give an 
approximation that is tolerably near. The theoretical possibility of repre¬ 
senting the distribution by a convergent series is not here considered, though 
the conditions of validity of the expansion are almost certainly satisfied. 
Hermite gives a method of obtaining the coefficients. Using the notation 

<£ (x) = g--a-(«* 3 + 2 hxy + by*)> ^ ^ _ e ~i- (bx^-2hxy+af)^ 

(d/dx) nl (d/dy) n (f) = 0 . U m , n , (d/dx) m (0/0?/) n yfr = yfrY mny 
U WjW and Y mn are both polynomials of degrees (m + n ) in x and y . If f(x, y) 
be expanded in the form n <j> U OTj n , it is proved that 

Mm, n 7=1 a rn, h| I f Cl/) "Vm, n^X dy y 

* ‘ Annali di Statistical 1883. 
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where u m%n is a numerical constant depending on m and n alone. Now it 
will he much more convenient in practice to obtain the coefficients by means 



of equating the successive 


moments 


corre- 


— 00 * —00 


sponding moments of the series. This is equivalent to Hermite’s process, 
Y m , n being a polynomial; so that n above is a linear function of the 
moments. 



ii r ^x m y n dx dy can be reduced by successive 


integrations by parts. If m< r, we have 

u r(t x m y n dx dy = (— ) m m\ 1 ly n u Ml} s dx dy = 0, 



whatever the values of s and and, similarly, the integral is zero if n<s, 
whatever the values of r and m. But if m = r, and n = s, 



.and if m > r and n > s } 

c rs x m y n dx dy = ( —) r+s (m !/(m—r) !) (n \j(n—s) !) x m ~ r y n ~ s dx dy, 




Using the notation cj> rs = j J (j>x r y s dx dy and fi r s = j j fx r y s dx dy, if we 
multiply the equation f(x , y) = 2 a rs , u rs by x r y* and integrate, we have 


fJL r $ — ^OO^rs““ ’^lO0r—l, s 4“ S$oi0r, s— l) 

+ {r (r— 1) a 2 o4>r-2 lS + r . ,-i + s (s— 1) &o20r,s-2}... , 


there being (r + 1) (s + 1) terms on the right, of which those are zero in 
which the sum of the subscripts of 0 is odd. 

Taking the first few equations, by choosing the origin at the mean of the 
distribution, and choosing the constants a, b , h so that 


/^20 • fill • fi02 * fiOO = 020 * 011 • 002 • 000, 


the constants a 0 1, «io, ^20, ^02 all become zero. By this means the normal 

surface that best fits the distribution is obtained, viz., fi oo0/0oo* 

But the fact that the contour lines in the example quoted are by no means 
concentric ellipses with centre at the mean shows that the fit is not yet 
sufficiently good. 

The next approximation to the distribution is 


/^000/000 + (^30^30 + 3/^21^21 4“ ^> fil 2 ^ h 2 + ^03^03) /6000- 
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The equations giving the coefficients of the fourth order terms are: 

/^40 = ^00^>40 + 4! ^40^00^31 = &00</>31 + 3 ! Ct3l<j)00> fi 22 = #00^22 +212! <X22<£()0» ©tC. 

Thus in terms of the moments fjb rs the coefficients are quickly obtained,, 
and Edgeworth's argument as to the rapidity of approximation applies 
equally well here if the standard deviations are sufficiently small. But the 
labour of obtaining the moments is very great, and a suitable example for 
testing the closeness of approximation is not available. What will be done 
here, therefore, is to show more or less generally how, in a simple case, contour 
lines can be obtained which are approximately of the type found in the 
example given above (see diagrams, Perozzo, loc . cit.). 

2. An Example of Contour Lines in Tivo Dimensions. —The example chosen 
here is that of <£ = e~^ x2 + and it will be assumed that there is symmetry 
about the line x = y. 

The most general first correction is a (%o + ^ 03 ) +/3 (^21 + ^ 12 ), where 
^30 = ( 3 ® — U21 = <f> y (1 — % 2 ), Ui 2 = cj)X (1 — y 2 ), it 0 3 = <£ (3 y — y s ). 

In particular, if /3 = 0, and a is negative, <*(^ 30 + % 3 ) vanishes on the 
straight line x + y = 0, and the ellipse x 2 —xy + y 2 = 3 ; on the positive 
side of x + y = 0 it is positive outside the ellipse, and negative inside, and 
vice versd on the negative side of x + y = 0. 

Thus the effect of the correction u (w 30 -f ^ 03 ), while not altering the mean 
or second moments, is to introduce skewness about x -f y = 0, the maximum 
ordinate being displaced in the negative direction along x = y. Outside the 
ellipse x 2 — xy + y 2 = 3 the surface is depressed on the negative side of 
x + y = 0 , and raised on the positive side. The zero-line of <j> + ot (uso + u os ) is, 
shown in the figure (curve A) for a = — 1/10. 

The contour lines on the positive part of the surface will now be oval curves, 
surrounding the point of maximum ordinate with their greatest width parallel 
to x + y, and crowding closely together in the neighbourhood of the zero- 
curve. These have not been shown in the figure to avoid complication. A 
variation in their form could be produced by introducing a term /3 (ui 2 + u 2 n) r 
but the effect of skewness is again the most noticeable. 

Consider now a fourth order correction, which for the end in view has been 
conveniently taken to be 

7 (^ 22 —w 3 i — W 13 ) = 7 {1— x 2 —y 2 +Gxy—xy (x 2 —xy + y 2 )}. 

The curve for which this is zero is drawn in the figure (curve B), and the 
portions of the plane are indicated in which it is positive and negative 
respectively. Thus the effect of this correction on the contour lines in the 
neighbourhood of the maximum ordinate is to contract them in the direction 
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.x — —y, and to elongate them in the direction x — y. A series of contour 
lines is drawn for 7 = 1 / 20 ; so that the lines shown are curves of equal 
frequency for the function 

{ 1 —(3 3 /3 ^ 3 + 0 3 /S^ 3 ) + yV (3 4 /cte 2 3y 2 —3 4 /3^3y 3 — 3 4 /3 ydx 3 )}e~~* + ^ /2) . 

The curves are distinctly similar in type to Perozzo’s curves, though no 
.attempt has been made at determining exactly what the coefficients should 
be. Probably a better approximation would be found with a rather larger 
coefficient for the third order correction, as the skewness is not so marked as 
it should be. 



Figure showing Contour Lines of the Surface 

z = {1 - X V (3 3 /oe 3 -1- 3 3 /Sy 3 ) + -.To (3 4 /3# 2 c?y 2 - d 4 /dx 3 dy - d 4 /dxdy 3 )} & 2+1 j 2 \ 

A— the curve ] 1 - T V (&/dx 3 + 3 3 dy 3 )} <*+*> = 0. 

B —the curve = {3 4 /3a; 2 3y 2 - 3 4 /3& 3 3y - 3 4 /3a% 3 } (2 ' 2+2/2) = 0. 
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Similar curves can clearly be obtained where there is correlation to be 
taken into account. For by making the change of variable 

7} = y), £-f ?; = /? (x + yj t e~^^ 2+2/2) 

is turned into a normal surface with correlation, while 8/3,r and 3/3y are 
linear functions of 3/3£ and 8/8/7, 


Eutectic Research . No. 1 ,—The Alloys of Lead and Tin . 

By Walter Bosenhain, B.A., B.C.E., with P. A. Tucker. 

(Communicated by Dr. B. T. Glazebrook, F.B.S. Beceived June 17,—Bead 

June 25, 1908.) 

(From the National Physical Laboratory.) 

(Abstract.) 

Attempts to prepare pure eutectic alloys of known constitution led to the 
discovery of discrepancies between the authors’ experiments and the data 
on lead-tin alloys published by Boberts-Austen. A complete redetermination 
of the equilibria of the lead-tin system was therefore undertaken, by both 
pyrometric and microscopical methods. Cooling-curves of the alloys taken by 
both inverse-rate and differential methods are given; these, together with the 
microscopic data, lead to the equilibrium diagram shown in the figure. This 
differs principally from that given by Boberts-Austen,* in that the eutectic 
point is placed at a concentration of 63 per cent, of tin instead of 69, that the 
eutectic line towards the lead end of the series terminates at a concentration 
close to 16 per cent, of tin, and that a series of transformations (along the 
line EFG) have been found in the solid alloys near the lead end of the series. 

The discrepancies of these results, as regards the eutectic composition, arise 
from the more delicate method employed in the present research. While alloys 
within 1 per cent, on either side of the true eutectic composition show no 
detectable difference in freezing- or melting-point, the presence of small 
amounts of excess of either constituent can be detected by the microscope, 
and in this manner the composition of the pure eutectic has been ascertained. 
As regards the solubility of tin in solid lead, it was found that the occurrence 
of the eutectic arrest-point in alloys containing less than 16 per cent, of tin 

* Boberts-Austen, ‘ Fifth Beport to the Alloys Besearch Committee of the Institutio 
of Mechanical Engineers, 5 1897. 









